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I. INTRODUCTION 

At a continuous transition only a few static exponents 
are enough to completely describe equilibrium critical be- 
haviour. When dynamic measurements are considered 
the critical behaviour becomes much richer, with further 
independent critical quantities having non-trivial critical 
exponents. 

It has been generally assumed that in a given dimen- 
sion all Ising Spin Glasses (ISGs) lie in a single univer- 
sality class. The Hamiltonian of the ISGs is given by 



(1) 



where Si = ±1 are the usual Ising spins and the nearest- 
neighbour couplings Jij are random variables. We have 
studied numerically critical dynamic behaviour for ISGs 
having near neighbour interactions in dimensions d = 3 
and 4 with Bimodal, Gaussian, or Laplacian distributions 
of the couplings. 1 We find that the dynamic exponents 
vary strongly and systematically from one distribution to 
another. One possible explanation for this could be that 
for ISGs the universality class depends on the form of the 
interaction distribution. 

In Section II we introduce the dynamical quantities, 
whereas in Section III we treat the important point of the 
reliable determination of critical temperatures. We then 
present our results for the spin glasses in four (Section 
IV) and three (Section V) space dimensions. Finally, 
Section VI contains our conclusions. 



II. DEFINITIONS 

In the field of dynamic critical measurements there is 
not complete consensus as to a universal convention for 
indicating protocols or exponents, particularly for the 
case of spin glasses. We will first define the convention 
that we will use, following 2 ^, and relate it to other stan- 
dard conventions. We will then summarize dynamic scal- 
ing properties. This section leans heavily on the review^. 



Throughout, we will implicitly consider only model A 
dynamics^ with single spin Glauber (or heat bath) up- 
dates in Ising spin glasses; the total time after quench is 
refered to as t. The basic dynamic protocol (there are po- 
tentially many others) consists in quenching the sample 
at time t = from T = oo to a temperature T, wait- 
ing (carrying out updates) for a time that will be called 
s; then at t = s and still at fixed T, either the phys- 
ical conditions are changed in some way (typically by 
switching on or switching off a small magnetic field h) or 
the instantaneous spin configuration at time s is simply 
registered to provide a reference state, without physical 
conditions being changed. There is then an observation 
period of further updates during which physical parame- 
ters are measured as functions of t. For convenience the 
time difference t — s is also denoted r. In alternative con- 
ventions the waiting time s is labeled t w , and r is labeled 
t. 

Among limiting conditions that can be profitably stud- 
ied for large samples are the condition s = (i.e. mea- 
surements start immediately on quenching), or alterna- 
tively s 3> t, a long waiting time condition after which 
the sample is in " quasi-equilibrium" . True equilibrium 
can be achieved if s is "long enough", a criterion that 
depends in a non-trivial way on the system, on the tem- 
perature T and on the sample size L. 

Many dynamic observables can be measured. Here we 
will concentrate on observations at criticality, T = T c , 
for the moment ignoring the question of how T c is to be 
estimated. 

A first fundamental definition is that of the dynamic 
critical exponent z c . At T c , the equilibrium autocorrela- 
tion relaxation time (with standard single spin updates) 
increases with sample size L as 

Tauto{L) ~ L Z < (2) 

where z c is the dynamical critical exponent. 

The two-time autocorrelation function is defined as 

N 



C(t ) 8) = ±[(£(S j (a)S j (t)) 



(3) 



where 



indicates the average over the thermal noise 
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and [■ • •] the average over the disorder (in alternative con- 
ventions C may also be written as q) . The critical scaling 
relation for C is 

C(t,s) = S - b f c (t/s). (4) 

In the quasi-equilibrium limit where s 3> t the critical 
scaling function f c (t/s) should follow the asymptotic be- 
haviour 

Ut/s) ~ [(t/s) lT b . (5) 

Thus in this limit we have 

C(t, s )~r- b (6) 

(in the alternative convention this is written as q(t) ~ t~ x 
with x = b) . For spin glasses the dynamic scaling relation 
governing b is£ 

b = (d - 2 + V )/2z c (7) 

where r\ is the static critical exponent. In the opposite 
limit when (t/s) — ► oo 

ut/s) ~ (v*r A ^ (8) 

or for s = 

f c (t)~t-^ z <. (9) 

A c /z c is related to the "initial slip" exponent 6 C (an in- 
dependent critical exponenli^) through 

9 C = d/z c - X c /z c . (10) 

The two-time linear autoresponse function is 

RQt, s) = [6(Si(t))/6h(s)} \ h=0 (t > s) (11) 

where h(s) is a time dependent conjugate magnetic field. 
The scaling equation for R is 

TR(t,s) = s-^f R (t/s) (12) 

with a — b for critical systems. The critical scaling func- 
tion f R (t/s) should follow the asymptotic behaviour 

fn(t/s) ~ (t/ s y x -/^ (is) 

when t/s — > oo. For short range inital correlations X R = 
A c . 

In simulations where the field h is applied at t = and 
switched off at t — s, the following integrated response is 
measured at times t > s: 

s 

p(t,s) = T J duR(t,u). (14) 
o 

This integrated response is directly related to the com- 
monly studied thermoremanent magnetization: 

M TR M(t,s) = hp(t,s). (15) 



Of further interest is the fluctuation-dissipation ratio 
X defined by 

X(t, s) = TR(t, s)/(SC(t, s)/Ss) = X{t/s). (16) 

In the quasi-equilibrium condition s>t the fluctuation- 
dissipation theorem holds and X — l. 4 When t/s — > 
oo, X takes a limiting value X^. For t s the ratio 
p(t, s) /C(t, s) also converges to this limit value. If the 
amplitudes A c and Ar arc defined in the same limit i>3 
by 

f c (t/s) = A c (t/s)- x </ z ° (17) 

and 

f R {t/s) - A R (t/s)- x ^^ (18) 

then2 

X oc = (A R /A c ){X c /z c -b)- 1 . (19) 

Finally, the dynamic spin glass susceptibility is measured 
through 

1 N 

Xne(t) = -[(Y / Sf(t)S^(t)) 2 } (20) 
J=l 

where a and /3 are two replicas of the same system relax- 
ing independently. The infinite time limit to the dynamic 
SG susceptibility is the equilibrium SG susceptibility for 
each size L, which at criticality increases with L as 

Xeq (L) ~ (21) 

where rj is the static critical exponent. The critical time 
dependence of the non-equilibrium spin-glass suscepti- 
bility for large samples after a quench to T c and with no 
anneal (s = 0) is& 

Xne(i) ~ i (2 -" )/2c = t h ' (22) 

where i is the time after quench and z c is again the dy- 
namical critical exponent. For convenience we have in- 
troduced an exponent h* — (2 — rj)/z c . 

Even for a canonical continuous transition such as that 
of the 2d Ising ferromagnet, where the static critical 
exponents are all known analytically and are rational 
numbers, the dynamic critical exponents can only be es- 
tablished numerically and have non-trivial values^*^ For 
Ising ferromagnets however, field theory (FT) epsilon ex- 
pansion estimates give reasonable agreement with numer- 
ical dynamic exponent estimates in dimensions 3 and 2. 4 

It is now well established that for standard systems 
that are in the same universality class not only the static 
exponents v, rj etc. but also the dynamic exponents z c , 9 C 
and Xrx, are all universal. The numerical data discussed 
in the following show that in each dimension for ISGs 
expected to lie in the same universality class the dynamic 
exponents vary strongly with the form of the interaction 
distribution. 
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III. ORDERING TEMPERATURES 

In order to obtain accurate and reliable simulation val- 
ues for the dynamic exponents an a priori requisite is to 
have reliable estimates for the ordering temperatures T c . 

High temperature series calculations give T c estimates 
which are not subject to finite size corrections and are 
thus intrinsically reliable, but whose accuracy is limited 
by the number of known terms in the series. Series es- 
timates can be extremely precise at high dimension but 
unfortunately they become progressively more inaccurate 
as the system dimension dropsi n i 12 ' 13 

There are a number of different ways in which to obtain 
estimates of ordering temperatures in ISGs through sim- 
ulations; for most of them one or more critical exponent 
estimates are also obtained simultaneously. Equilibrium 
simulations can provide accurate data but necessarily on 
samples which are of small or moderate size L, and the 
reliability and precision of the T c estimates are finally 
limited by the need to extrapolate to large L to elimi- 
nate corrections to finite size scaling whose importance 
depends on the system being studied, on the parameter 
being measured, and on the maximum range of sample 
sizes that can be equilibrated with the computing facili- 
ties available. Corrections to scaling are subtle even for 
the canonical Ising ferromagnet:si2ii£ii£i where the lead- 
ing corrections include both "irrelevant operator" and 
"analytic" contributions, while for ISGs basic guide-lines 
are lacking so that one must rely on empirical observa- 
tions. In principle it should be possible to use the onset 
of deviations from strict critical behaviour to monitor 
T c ; for instance ln(x(X, T)) against ln(L) curves bend 
downwards/upwards for T greater/less than T c . In the 
SG context this approach has been rarely used as the 
upbending below T c is weak, at least in 3d. The equi- 
librium finite size scaling simulation techniques, which 
are most often relied on to estimate T c in ISGs, are 
based on measurements of parameters that are dimcn- 
sionless and take on an L-independent value at T c for 
large L. Well known examples are the Binder moment 
ratio^i g(L, T) and the second moment correlation length 
ration T) /L. Plots of g(L, T) or T) /L as func- 
tions of T for fixed L have a unique " crossing point" at 
a temperature which is equal to T c in the limit of large 
L. These methods require strict thermal equilibration at 
each size L; also the exact position of the large L limit 
crossing point may be masked up to quite large L by cor- 
rections to scaling (see for instance^ for the case of the 
canonical Ising ferromagnet in three dimensions). 

Simulations become heavier with increasing d simply 
because for given L the number of spins is L d , but this 
effect is compensated by the fact that z c tends to drop 
with increasing d. Furthermore crossing points become 
better denned at higher d, and it turns out that correc- 
tions to finite size scaling become weaker as d increases. 
On balance it is in fact easier to estimate T c reliably by 
equilibrium simulations at d = 4 (and above) than at 
d = 3. 



An alternative simulation technique which we will rely 
on below is to combine static and dynamic measurements 
to estimate T c by consistency^ This has the advantage 
of using two dynamic measurements which do not re- 
quire equilibration and which have negligible corrections 
to scaling, together with equilibrium spin glass suscepti- 
bility measurements which do require equilibration but 
which are less sensitive to corrections to scaling than are 
Binder moment ratio or correlation length ratio measure- 
ments. A range of putative T c values T* are chosen, and 
three measurements are made at each T* : 

- the effective dynamic exponent b{T*) from large 
L quasi-equilibrium measurements with s S> r using 
Eq. ®, 

- the effective dynamic exponent h*(T*) from large L 
measurements of the dynamic SG susceptibility Eq. 122|l . 
and 

- the effective static exponent r/(T*) from equilib- 
rium SG susceptibility finite size scaling measurements, 
Eq. gl}. 

As these three parameters are linked at T c through the 
two exponents r\ and z c , Eq. Q, Eq. (2} and Eq. (22J, 
there is a consistency condition which holds at and only 
at T* — T c . There are different ways to implement this 
condition. We can first use the equilibrium and dynamic 
SG susceptibility results together to obtain a set of val- 
ues of z(T*) at each T*,— These values are to good pre- 
cision independent of corrections to scaling, and the set 
of z(T*) extends from above to below T c . Secondly, from 
the effective r)(T*) and b(T*) one can derive a second ef- 
fective z, z+(T*) = (d- 2 + r](T*))/2b{T*). At T = T c 
the consistency condition is simply z(T*) — z + (T*). The 
measured values of the parameters r/(T*) and z(T*) at 
this unique temperature correspond to the true critical 
exponents r\ and z c . This method is rather insensitive to 
corrections to finite size scaling. 

The three ISG distributions which will be considered 
here are the random Bimodal, Gaussian, and Laplacian 
near neighbour interaction distributions on (hyper)cubic 
lattices. The explicit normalized distributions are 

P B (Jij) = [S(Jij -J)+ S(J tJ + J)]/2, (23) 

P G (J l3 ) = exp(-4/2J 2 )/(jV2^) (24) 

and 

P L {J l3 ) = exp(-\/2 | J l3 /J |)/(J\/2) (25) 

respectively. The distributions are symmetric about zero 
and are normalized in such a way that < Jfj > / J 2 = 1 . 

IV. DIMENSION 4 

We will first consider explicitly dimension 4. T c val- 
ues for a range of interaction distributions including the 
three cases that concern us here were obtained by simu- 
lations using the consistency methodi For the Bimodal 
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FIG. 1: Scaling behaviour of the autocorrelation function 
C(t,s) with a Laplacian distribution of the couplings in (a) 
four and (b) three dimensions. The best data collapses yield 
the values of b given in Tables HI and HT1 



and Gaussian distributions the values were fully consis- 
tent with and were as accurate or more accurate than 
other simulation estimates using alternative simulation 
techniques ^2i2ii22i2&i2ii No other result appears to have 
been reported for the Laplacian distribution. Essentially 
negligible corrections to scaling can be seen in the data 
for any of the simulation techniques at this dimension; 
for instance, the Binder parameter crossing points are 
well defined and appear to be independent of L to within 
high numerical precision. High temperature series esti- 
mates for the Bimodal caseii*i2ii^ and for other interac- 
tion distributions^ are in excellent agreement with the 
simulation estimates; an overview of the data is given 
iriSi. The overall agreement between the complementary 
approaches means that in d = 4 the consistency sim- 
ulation technique^ is validated. Hence the T c values, 
together with the associated 77 and z c critical exponent 
values from the consistency method, can be taken as re- 
liable. 

Non-equilibrium measurements of the two-time auto- 
correlation function and the two-time linear autoresponsc 
function were made at the temperatures corresponding 
to the T c values estimated from the consistency method. 
Large systems containing 20 4 spins were simulated us- 
ing the standard heat-bath algorithm. The systems were 
prepared initially in a completely disordered state and 
then quenched down to T c at time t = 0. For the compu- 
tation of the thermoremanent magnetization an external 
field with strength h = 0.05 was applied between t — 
and t = s with s varying from 25 to 400. 

Figure^ and Figure [2] summarize our findings for the 
four-dimensional systems. The expected dynamical scal- 
ing behaviour of the autocorrelation function is illus- 
trated in Figure for the case of a Laplacian distribu- 



FIG. 2: Out-of-equilibrium quantities in dimension four, (a) 
Time depencence of the autocorrelator with s = 0. The three 
different distributions (b: Bimodal, g: Gaussian, 1: Lapla- 
cian) yield different exponents for the power-law behaviour 
observed at long times, (b) Fluctuation-dissipation ratio X 
as function of s/t. In the limit s/t — > one obtains the limit 
value Xao with different values for the different distributions. 



tion of the couplings. Plotting C(t, s) as a function of 
t/s for various values of the waiting time s, an excellent 
data collapse is achieved for the value b = 0.140(3). De- 
viations from this scaling behaviour are only obvious in 
the regime t — s < s, i.e. outside of the dynamical scaling 
regime. A similar good data collapse^ is obtained for the 
other distributions, see Tableland Refi£. It is worth not- 
ing that the values of b obtained in these non-equilibrium 
simulations agree with those derived from the quantities 
77 and z c via Eq. Q. Equilibrium and non-equilibrium 
simulations therefore consistently yield for ISGs critical 
quantities depending on the form of the distribution of 
the couplings. 

In Figure [3 we discuss truly non-equilibrium quanti- 
ties which can not be expressed solely by equilibrium 
quantities. As shown in Figure^ plotting lnC(i,0) ver- 
sus ln/j results in straight lines in the long time limit, 
in agreement with the expected power-law behaviour 
(PI). The slopes of these lines yield the exponent \ c /z c . 
Again, this quantity, supposed to be universal, shows a 
clear dependence on the chosen distribution, see Table 
H] Finally, Figure Eb displays the temporal evolution of 
the fluctuation-dissipation ratio Hlfi|l which in the limit 
t/s — ► 00 yields the limit value X 00l again supposed to 
be universal. It is obvious from this plot that the value 
of is different for the three distributions considered 
in this work. 

Values for the various parameters corresponding to the 
present three distributions are shown in Table [I] The 
amplitude ratio Ar/A c has been derived from Eq. 119|l . 

By inspection of the results in Table I it can be seen 
that the equilibrium critical exponent r\ together with 
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TABLE I: 


Parameter estimates in 


dimension 4. 


parameter 
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all the dynamic critical exponents vary strongly from 
one distribution to another. Apparent non-universality 
of critical exponents obtained from simulation data has 
in the past been ascribed to a consequence of errors 
in the estimation of critical temperatures or to a lack 
of care in allowing for corrections to finite size scal- 
ing. In the present case, the values of the ordering 
temperatures^ have been validated by high temperature 
series calculations^ and internal evidence, inherent to 
the consistency method described in Section III, shows 
that corrections to finite size scaling are negligible in this 
dimension. In addition, the dynamic parameters are ob- 
tained from simulations on large samples and can be con- 
sidered virtually free of finite-size corrections to scaling. 
The values of each of the dynamical critical parameters 
are insensitive to the precise value of the ordering tem- 
perature so even if there were small errors in the assumed 
values of the ordering temperatures the effects on the dy- 
namic critical parameter estimates would be negligible. 

At this point some remarks on the reliability of the 
extraction of critical exponents from out-of-equilibrium 
simulations are in order. The whole approach is based on 
the assumption that, in the time range in which the ex- 
ponents are determined, the dynamical correlation length 
£(t) increases as a simple power-law 



(26) 



where z c is the dynamical critical exponent. For critical 
ferromagnets this growth law and the resulting dynamical 
scaling set in already after a few time steps. However, as 
spin glasses are characterized by a large value of z c one 
may wonder whether this simple growth low prevails for 
the times we have accessed in our simulations or whether 
a more general growth law of the form 



£(t) = at 1 ' 2 " + bt 1 ' 2 



(27) 



with a sizeable finite-time correction is observed. 

In fact, the growth of the dynamical correlation 
length at the critical point of different three- and four- 
dimensional spin glasses has been intensively investigated 
in the recent past. We mention here the Ising spin glass 
with a Bimodal 29,30 or a Gaussian 21,30 distribution of the 
couplings, the gauge glass with a Gaussian distribution 21 , 
the XY spin glass with a Bimodal distribution^ or the 



FIG. 3: Dynamical exponents z (filled circles) and z + (open 
squares) as function of trial temperatures T* for ISGs in di- 
mension three with (a) Bimodal, (b) Gaussian and (c) Lapla- 
cian distribution of the couplings. 



Heisenberg spin glass with a Gaussian distribution 3 ^. All 
these studies reveal that for times t > 20 the increase of 
the dynamical correlation length in various spin glasses 
(including some of the cases we consider in this work) 
is given by the simple power-law H2(jfl . No finite-time 
corrections of the form (|27ll have been observed. Addi- 
tional support for the growth law comes from the 
perfect dynamical scaling behaviour of two-time quanti- 
ties (as for example the autocorrelation function shown 
in Figure [Q, as a sizeable finite-time correction would 
completely spoil the observed data collapse. 



V. DIMENSION 3 

In dimension d = 3 the overall situation is rather 
less satisfactory; high temperature series valuesii*^ be- 
come imprecise for the Bimodal case and none have been 
reported for the other distributions. The Binder ratio 
method becomes delicate because the gh{T) curves only 
fan out weakly at low temperatures making the limit- 
ing crossing point difficult to identify and very sensi- 
tive to corrections to scalingiii The correlation length 
ratio appears to suffer from strong corrections to scal- 
ing, especially at low LiS For the Gaussian distribu- 
tion there is a general consensus as to the value of T c 
from different simulation estimates For the Bi- 

modal distribution, published T c estimates are much 
more scatterec&i^iiSi 3 ^! which we ascribe to difficulties 
related to corrections to finite size scaling. For the Lapla- 
cian distribution we are not aware of other published es- 
timates. 

We will rely on the value from the consistency method 
because, for the reasons outlined above, this technique 
is much less sensitive to problems of corrections and 
because it has given excellent agreement with the high 
temperature series values in d = 4. Data from the 
consistency method are presented in Figure 3. With 
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TABLE II: Parameter estimates in dimension 3. 



parameter 


Bimodal 


Gaussian 


Laplacian 


T c 


1.19(1) 


0.92(1) 


0.72(2) 


Zc 


5.7(2) 


6.2(1) 


8.6(2) 


V 


-0.22(2) 


-0.42(2) 


-0.55(2) 


b 


0.056(3) [2] 


0.043(1) [2] 


0.032(2) 


Xc/Zc 


0.362(5) [2] 


0.320(5) [2] 


0.259(2) 




0.12(1) [2] 


0.09(1) [2] 


0.055(2) 


6c 


0.165(25) 


0.165(13) 


0.090(10) 


A R /A C 


0.037(4) 


0.025(4) 


0.0125(7) 



0.20 




ln(t) s/t 

FIG. 4: The same as Figure 2, but now in dimension three. 
Whereas the insert in (a) shows the time dependence of the 
effective exponent (A c /z c )e// 1281 . the insert in (b) displays 
the long time behaviour of X. 

sets of trial temperatures T* we plot for each system 
z{T*) and z+{T*) against b(T*). z(T*) is derived from 
a comparison of the equilibrium and dynamic SG sus- 
ceptibility results^ at each T*, and z+'(T*) = (d - 2 + 
rj(T*))/2b(T*). b(T*) is measured from the autocorrela- 
tion function decay in quasi-equilibrium as defined above, 
and the effective rj(T*) is obtained from equilibrium finite 
size SG susceptibility measurements. At T c consistency 
of the various exponents dictates that z(T*) = z + {T*). 
The values of T c together with the exponents z c and rj ob- 
tained are given in Table II (the values are more precise 
than those given using the same method inSi because of 
improved equilibrium susceptibility data 33-38 ). It is in- 
teresting to note that Migdal-Kadanoff estimates of T c 
for the different distributions with d = 3 and the MK 
parameter b = 222a£ are strikingly similar to the values 
given here. 

For the non-equilibrium simulations in three dimen- 
sions we considered systems with 50 3 spins and wait- 
ing times s < 1600. The expected scaling behaviour of 
the two-time quantities is again observed, see Figure^) 
for the autocorrelation of the Laplacian distribution. As 
shown in Figure 0] (see Table [nj one observes that the 
values of \ c / z c and Xoo also depend in three dimensions 
on the form of the distribution function of the couplings. 
The insert in Figure 0Ji displays the effective exponent 

(X c /z c ) eff = -(ln(C(20t,0)-ln(C(i,0))/(ln(20t)-ln(t)) 

(28) 

as a function of 1/i. In all cases this effective exponent 
rapidly reaches a constant, distribution dependent, value. 

Once again the estimates of the various static and dy- 
namic critical exponents vary considerably from distribu- 
tion to distribution. The sense of the variations is system- 
atic and is the same as in dimension 4 : with increasing 
kurtosis of the distribution, T c drops, z c increases, rj be- 



comes more negative, and the dynamic exponents either 
expressed as X c /z c and X^ or as 9 C and Ar/A c all drop. 



VI. DISCUSSION 

For standard continuous transitions the renormaliza- 
tion group theory provides a comprehensive explanation 
of critical behaviour and in particular of the strict iden- 
tity of exponents of all systems within each universality 
class, the class being defined by a restricted list of pa- 
rameters which includes the physical dimension d and the 
number of order parameter components n. The univer- 
sality covers not only equilibrium exponents but extends 
to the whole family of dynamic exponents. This uni- 
versality reflects the fundamental principle that within 
each class, the details of the physics at the local level do 
not affect the large scale behaviour which determines the 
critical exponents. 

It has been widely assumed that in a given dimension 
all ISGs fall in the same universality class. However it 
should be noted that the critical behaviour of spin glasses 
is qualitatively very different from that of standard sys- 
tems such as ferromagnets. The upper critical dimension 
is 6 rather than 4, and below the upper critical dimension 
the specific heat exponents are strongly negative so there 
is no specific heat peak or cusp. Field theory (which pro- 
vides the well known e expansion development at stan- 
dard continuous transitions) has proved intractable in the 
ISG context below the upper critical dimension d = 641 
Already at d = 5 and d = 4, numerical values of the equi- 
librium exponents obtained from summing the known 
leading terms to order three in the ISG e expansion^ are 
very different from estimates using the high temperature 
series methodiSii^ or simulations. For the Ising ferromag- 
net at d = 3 and so e = 1, the FT development in e to 
third order is accurate to better than 0.001. 43 This is in 
total contrast to the situation for the Bimodal ISG where 
at d = 5 (so again at e = 1) the FT sum to third order in e 
gives r](d = 5) = 1.6897, strikingly different from the high 
temperature series valued ri(d = 5) = —0.38(7) and the 
simulation estimate T](d = 5) = — 0.39(2).! This implies 
that a sum including many further terms (oscillating in 
sign) would be needed to finally obtain stable and accu- 
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rate FT predictions. In practice, establishing such a sum 
seems entirely ruled out, but the question remains open 
as to whether the necessary quasi-cancellations among 
the unknown higher order FT terms could depend on pa- 
rameters such as the lattice structure or the form of the 
interactions. 

In our simulations of the Ising spin glasses we have 
found that, in most cases, the differences between the 
exponent estimates for the different systems are much 
larger than the statistical error bars. This is especially 
obvious for the amplitude ratio Ar/A c which is supposed 
to yield clearly distinct values for different universality 
classes, similar to what is observed for static ampltudc 
ratios at equilibrium. Extrapolations in each dimension 
suggest that if the interaction distribution is modified 
and T c decreases, all the exponents studied vary system- 
atically in such a way that z c increases strongly, rj tends 
towards a value near 2 — d (which is the strict limiting 
value for 77 in each dimension when T c — 0) while 
and 9 C tend to near zero. The data as they stand are thus 
compatible with exponents each varying continuously to- 
wards a T c — limit as T c is driven lower by a widening 
of the interaction distribution (increasing kurtosis). 

Our non-equilibrium simulations can in principle be 
subjected to only three different sources of systematic 
errors: (1) the values of the critical temperatures are 
erroneous, (2) the sizes of the samples are too small, or 
(3) the time range of our runs is insufficient. Let us 
address these three different points. 

The values of the critical temperatures we use have 
been estimated with a technique combining equilibrium 
and non-equilibrium measurements^ and are in excellent 
agreement with independent simulation estimates. In ad- 
dition, T c in the four-dimensional systems have been con- 
firmed recently^ by high temperature series estimates. 
This agreement validates the technique used for the eval- 
uation of T c and also, as the method relies on a consis- 
tency argument, indirectly the values of z c and rj obtained 
in the same simulations. 

The data we have discussed in this paper concern prin- 
cipally dynamic exponents in ISGs. These measurements 
have the advantage of not requiring strict equilibration, 
which in turn permits studies on samples so large that 
they cannot be conveniently equilibrated at or near crit- 
icality. As the sample sizes L are much larger than the 
maximum of the correlation length attained dur- 
ing the simulations, the measurements are always taken 
in the infinite sample size limit and are not hampered 
by finite-size effects. We have thus been able to min- 



imize one of the major sources of systematic errors in 
numerical simulations, namely corrections to finite size 
scaling. This is very similar to what is observed when 
studying critical ferromagnets, as for example the three- 
dimensional Ising model 44 , where no notable finite-size 
corrections to scaling are encountered in non-equilibrium 
simulations, in contrast to equilibrium simulations where 
these corrections to scaling can be very strong. 

Whereas finite-size effects are well controlled in non- 
equilibrium simulations of large systems, this is not im- 
mediately obvious for finite-time corrections. Indeed, for 
spin glasses the critical dynamical exponent z c takes on 
very large values, wheh might raise some doubts whether 
the simple power-law increase (|26[1 of the dynamical cor- 
relation length is valid in the time range we accessed. 
However, various investigations of spin glasses in differ- 
ent dimensions and with different distributions of the 
couplings in the recent past have found that in general 
the increase of is completely described by (12(it for 
times t > 20. We can therefore be confident that the 
dynamical scaling approach, underlying our estimates of 
the critical quantities, is also valid for spin glasses and 
that the values of the equilibrium and non-equilibrium 
critical quantities we obtain are reliable. Furthermore, 
as shown in Figure 01 the run times of our simulations 
are clearly sufficient for the dynamic parameters to take 
up their limiting values. 

Having now excluded the most probable sources of sys- 
tematic errors we interpret our numerical data as strong 
evidence that in spin glasses critical quantities do depend 
on the exact form of the distribution of the couplings. 
Of course, as we only provide numerical evidence, we 
can not completely exclude that corrections coming from 
other sources could have some impact on the values of 
the critical exponents. 
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